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We investigate how stabilizer theory can be used for constructing sufficient conditions for entan- 
glement. First, we show how entanglement witnesses can be derived for a given state, provided 
some stabilizing operators of the state are known. These witnesses require only a small effort for 
an experimental implementation and are robust against noise. Second, we demonstrate that also 
nonlinear criteria based on uncertainty relations can be derived from stabilizing operators. These 
criteria can sometimes improve the witnesses by adding nonlinear correction terms. All our criteria 
detect states close to Greenberger-Horne-Zeilinger states, cluster and graph states. We show that 
similar ideas can be used to derive entanglement conditions for states which do not fit the stabilizer 
formalism, such as the three-qubit W state. We also discuss connections between the witnesses and 
some Bell inequalities. 
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I. INTRODUCTION 

Entanglement lies at the heart of quantum mechan- 
ics and plays also a crucial role in quantum information 
theory. While the properties of bipartite entanglement 
are still not fully understood, the situation for multi- 
partite entanglement is even more unclear, since in the 
multipartite setting several inequivalent types of entan- 
glement occur. However, entangled states of many qubits 
are needed for quantum information tasks such as mea- 
surement based quantum computation error correc- 
tion 0, Q or quantum cryptography Q , to mention only 
a few. Thus it is important both theoretically and ex- 
perimentally to study multipartite entanglement and to 
provide efficient methods to verify that in a given exper- 
iment entanglement is really present. 

In this paper we will apply the stabilizer theory 0, Q 
for entanglement detection. This theory already plays a 
determining role in quantum information science. Its key 
idea is describing the quantum state by its so-called stabi- 
lizing operators rather than the state vector. This works 
as follows: An observable Sk is a stabilizing operator of 
an iV-qubit state \t/j) if the state \ip) is an eigenstate of 
Sk with eigenvalue 1 

«> = l# (1) 

Many highly entangled iV-qubit states can be uniquely 
defined by N stabilizing operators which are locally mea- 
surable, i.e., they are products of Pauli matrices. 

The main result of the present paper can be formu- 
lated as follows: If one has a given state \tp) and has 
identified some of its stabilizing operators, then it is easy 
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to derive entanglement conditions which detect states in 
the proximity of \ip). So looking for stabilizing opera- 
tors should be the first step in order to detect entangle- 
ment. All the conditions presented are easy to imple- 
ment in experiments and are robust against noise. We 
mainly consider criteria detecting entanglement close to 
Greenberger-Horne-Zeilinger (GHZ) cluster Q and 
graph states • We use different types of entanglement 
conditions, they may be linear criteria, such as entan- 
glement witnesses, or nonlinear criteria, based on uncer- 
tainty relations. In this way we complete our results of 
Ref. [3 , where witnesses for the detection of multipartite 
entanglement in the vicinity of GHZ and cluster states 
were derived. Note that the stabilizer formalism also al- 
lows to derive Bell inequalities, this has been recently 
investigated in Ref. |9(. We do not aim to derive Bell 
inequalities here, since they are used to rule out local 
hidden variable models, a notion independent of quan- 
tum physics. However, some of our constructions exhibit 
close connections to Bell inequalities, and this will also 
be discussed. 

Our paper is organized as follows. Since GHZ states 
are the most studied stabilizer states, we use mainly them 
to explain our ideas, the generalization to other stabilizer 
states is then usually straightforward. So we start in Sec- 
tion II by recalling the basic facts about the stabilizing 
operators of GHZ states. Then we present a method for 
obtaining a family of entanglement witnesses for detect- 
ing entanglement close to GHZ states. First we present 
witnesses detecting any (i.e., even partial or biseparable) 
entanglement. Then we present witnesses which detect 
only genuine multi-qubit entanglement. We discuss some 
interesting connections to Bell inequalities. In Section 
III we present witnesses for cluster and graph states. We 
also consider detecting entanglement close to given mixed 
states. Finally, we present entanglement witnesses for a 
W state. It is of interest since the W state does not fit the 
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stabilizer framework. We show that our ideas can still be 
generalized for this case. In Section IV we present non- 
linear entanglement conditions in the form of variance 
based uncertainty relations. It turns out that they can 
often improve the witnesses by adding nonlinear terms. 
In Section V we present entanglement conditions which 
are based on entropic uncertainty relations. Finally, in 
the Appendix we collect some basic facts about the sta- 
bilizer formalism and present some technical calculations 
in detail. 



state is fully separable, i.e., it can be written as 



= ^PiQi ® Qi ® ■■■ <8> Q. 



AN) 



(5) 



with pi > 0,J2iPi = !• 

Before presenting entanglement witness operators, let 
us shortly recall their definition. An entanglement wit- 
ness W is an observable which has a positive or zero 
expectation value for all separable states, and a negative 
one on some entangled states [Toj| : 



II. GHZ STATES AS EXAMPLES OF 
STABILIZER STATES 

Let us start by introducing the stabilizer formalism 
using the example of GHZ states. An A-qubit GHZ state 
is given by 



\GHZ N ) = -L(|0) 



II 



(2) 



Besides this explicit definition one may define the GHZ 
state also in the following way: Let us look at the ob- 
servables 



S\ 



(GHZ N ) 



N 
k=l 



s 



(GHZ N 



) ,_ Z {k-i) z {k) for fe = 2, 3, AT, (3) 



where X^ k \ Y^ k \ and denote the Pauli matrices act- 
ing on the fc-th qubit. Now we can define the GHZ state 
as the state \GHZn) which fulfills 



S 



(GHZ N ) 



GHZ N ) = \GHZ 



(4) 



for k = 1,2,..., A. One can straightforwardly calculate 
that these definitions are equivalent and the GHZ state 
is uniquely defined by the Eqs. 10}. From a physical point 
of view the definition via Eqs. |@} stresses that the GHZ 
state is uniquely determined by the fact that it exhibits 
perfect correlations for the observables Si HZ . 

Note that \GHZ N ) is stabilized not only by s[ GHZn \ 
but also by their products. These operators, all hav- 
ing perfect correlations for a GHZ state, form a group 
called stabilizer 0. This 2 Ar -element group of operators 



will be denoted by S^ GHZn \ SY 



(GHZ N ) 



are the genera- 



tors of this group which we will denote as S^ GHZn ^ — 



{GHZ h 



i{GH Z]\ 



, Q (GHZ N ) 



). For more details on 



the stabilizer please see Appendix A. 



Tr{Wg) 



> for all separable states g s . 
< for some entangled states g e 



(6) 



Thus a negative expectation value in an experiment sig- 
nals the presence of entanglement. 

In this paper we will construct a family of entangle- 
ment witnesses using the elements of the stabilizer. Wc 
will call these stabilizer witnesses. 



B. Ruling out full separability 

In the following stabilizer witnesses will be used to de- 
tect entanglement close to GHZ states. We will construct 
witnesses of the form 

W := c t - S { k GHZN) -S\ GHZn \ (7) 
where S^^ Zn ^ are elements of the stabilizer group, 



'k/l 



cq := max \{Sj 



qev 



(GHZn) qi^GHZn)^ 



(8) 



and V denotes the set of product states. Since the set 
of separable states is convex, cq is also the maximum for 
mixed separable states of the form Eq. JSJ). 

Clearly, if we want to detect entangled states with W, 
we have to choose g l Y GHZN '> anc j gi GHZ ff) the 



of (S k 



(GHZ h 



S^ GHZn ^) for entangled quan- 



maximum 

turn states is larger than the maximum for separable 
states. Whether this condition holds depends on the 



and S, 



(GHZ N ) 



commute 



question whether the S k 
locally: 

Definition 1. Two correlation operators of the form 



K = A (1) <g> K {2) 
L = (g)L (2) € 



...® A-W 



(9) 



commute locally if 
for every n G {1,2, A} : K^L^ = L^K (n \ (10) 



A. Witnesses for stabilizer states 

In order to show that a given state contains some en- 
tanglement, we have to exclude the possibility that the 



Using Definition 1, we can make the following statement: 
Observation 1. Two multi-qubit correlation operators, 
K and L, commute locally iff there is a pure product state 
among their common eigenstates. 

Proof. K and L commute locally iff for all n € 
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{1,2, ...,N} there are two vectors, \<f) n ) and \(j>n), which 
are common eigenstates of and L^ n \ Thus \ip) — 

|0i) ® |02) ® ■•■ ® |0at) is a common eigenstate of K and 
L. □ 
Hence it follows that if S < ^ HZn * > and S^ 0112 "^ commute 



(GHZ N ) 



s. 



(GHZf, 
I 



for 



locally then the maximum of (5 ( 
separable and entangled states coincide. 

After these considerations, we construct our witness 
from two locally non-commuting stabilizing operators: 
Theorem 1. A witness detecting entanglement around 
an N-qubit GHZ state is 



yy(GHZ N ) ,_ t _ S (GHZ N ) _ S (GHZ N )^ 



(11) 



where m = 2,3, N. 

Proof. The proof is based on the Cauchy-Schwarz in- 
equality. Using this and (X^) 2 + (Z^) 2 < 1, for pure 
product states we obtain 

{S (GHZ N)) + {S {GHZ N)) = 

= (lW)(l( 2 ))...(lW) + (Z^ m -^)(Z^) 

< nx^-^l ■ \(x {rn) )\ + {(z^-^i ■ \(z {m} )\ 

< y^™- 1 )) 2 + (z(™- 1 )) 2 y / (x(™)) 2 + (z(™)) 2 

< I- (12) 

It is easy to see that the bound is also valid for mixed sep- 
arable states. This proof can straightforwardly be gener- 
alized for arbitrary two locally non-commuting elements 
of the stabilizer. □ 
Witnesses can be constructed with more than two ele- 
ments of the stabilizer as 



GHZjs 



_ -,. q{GHZnt) Q (GHZ N ) Q {GHZ N ) Q (GHZ N ) 
— JL D 1 D m D 1 D m 

(13) 

for m = 2,3, N rule out full separability. This can be 
proved in a similar calculation as in Theorem 1, using 
the fact that (X^) 2 + (F«) 2 + (Z«) 2 < 1. Later, in 
Sec. IV we will see how these conditions can be improved 
by adding refinement terms which are quadratic in the 
expectation values. 



C. Criteria for witnesses 

Having derived the first entanglement witnesses, it is 
now time to ask, whether they are really useful witnesses. 
In fact, there are, for a given state, always infinitely many 
witnesses allowing the detection of this state. However, 
some of them are more useful than others. Two crite- 
ria for the usefulness of witnesses are of interest: Firstly, 
it is important to characterize how much a witness tol- 
erates noise. Secondly, it is crucial to know, how much 
experimental effort is needed to measure the witness. 

Witnesses are usually designed to detect entangled 
states close to a given pure state \*S?). From a practi- 
cal point of view it is very important to know, how large 



neighborhood of |\E r ) is detected as entangled. This can 
be characterized in the following way: 
Definition 2. Let the density matrix of the state ob- 
tained after mixing with white noise be given by 

f?(p„oisc) ^Pnoisc^T + (l-Pnoise)|*)(#|- (14) 

Here p n oise determines the ratio of white noise in the mix- 
ture. Then, the robustness to noise for a witness W is 
determined by the maximal noise ratio for which it still 
detects £>(p n oise) as entangled. 

It is easy to see that witness W detects £>(p n oise) as 
entangled if p no isc < Piimit where 



Plimit : = 



-(GHZ N \W\GHZ 



2- N Tr{W) - (GHZ N \W\GHZ_ 



(15) 



To give an example, for witnesses of the form Eq. J7J we 
have piinut = (7 - l)/7 where 



{GHZ N \sf HZN) 



S. 



(GHZ N ) 

I 



GHZ 



N 



maxgg-p [(S 



(GHZ N ) 
k 



S, 



{GHZ* 



(16) 



Clearly, the maximum of the numerator is 2, while the 
minimum of the denominator is 1. Thus the maxi- 
mum noise tolerance achievable by a witness of the form 
Eq. is given by p limit = 1/2. Witnesses W { ^ HZn) 
from Eq. (|ll|l have exactly this noise tolerance thus they 
are optimal among stabilizer witnesses with two corre- 
lation terms. Witnesses Wm HZN ^ from Eq. (|13fl detect 
£>(Pnoise) as entangled if p no ise < 2/3, thus, they are more 
robust against noise. It can be proved that Eq. (|13|) is op- 
timal among stabilizer witnesses with three correlations 
terms from the point of view of noise tolerance. 

The experimental effort needed for measuring a witness 
can be characterized by the number of local measurement 
settings needed for its implementation [Tlj : 
Definition 3. The local decomposition of a witness is 
defined as follows: Any witness can be decomposed into 
a sum W = Y)^ M.j where each of the terms Mi can 
be measured by a local measurement setting. One local 
setting C = {0^}fcLi consists of performing simultane- 
ously the von Neumann measurements on the corre- 
sponding parties. By repeating the measurements many 
times one can determine the probabilities of the 2 N dif- 
ferent outcomes. Given these probabilities it is possible to 
compute all two-point correlations (O^OW), three-point 
correlations (O^O^O^), etc. 

It is easy to see that measuring only one setting 
is not enough for entanglement detection. For mea- 
suring all the witnesses Wm HZn ^ given Eq. 1111) for 
m = 2,3, ...,N, two measurement settings are required 
for an implementation, namely {X^\ X^ 2 \ X^ N '} 

and {ZW,Z®,...,ZW}. The witnesses w£ BZlf) 
from Eq. (|13[) require a measurement of three set- 
tings: {X^\X( 2 \...,X^}. {FW,^ 2 ), and 
{ZV,ZW,...,ZW}. 
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D. Detecting genuine multipartite entanglement 

Up to now we considered witnesses which detect any 
(even partial) entanglement. However, for multipartite 
systems there are several classes of entanglement. The 
most interesting class of entangled states are the genuine 
multipartite entangled states. These are defined as fol- 
lows. Let us assume that a pure state \ip) on an TV qubit 
system is given. If we can find a partition of the N qubits 
into two groups A and B such that is a product state 
with respect to this partition, 

|V>> = |^|X>B (17) 

then we call the state \tp) biseparable (with respect to 
the given partition). Note that the states \(/))a and \x)b 
may be entangled, thus the state \ip) is not necessarily 
fully separable. According to the usual definition a mixed 
state is called biseparable iff it can be written as a convex 
combination of pure biseparable states 

= 5>#;)<^| (18) 

i 

where the l^) may be biseparable with respect to differ- 
ent partitions. If a state is not biseparable then it is called 
genuine multipartite entangled. In experiments dealing 
with the generation of entanglement in multi-qubit sys- 
tems it is necessary to generate and verify genuine mul- 
tipartite entanglement, since the simple statement "The 
state is entangled" would still allow that only two of the 
qubits are entangled while the rest is in a product state. 

Witnesses for the detection of genuine multipartite en- 
tanglement have already been used experimentally for en- 
tanglement detection close to three-qubit W states and 
four-qubit singlet states 0. They used the projector on 
the state to be detected as an observable to detect en- 
tanglement. For GHZ states, the projector-based witness 
reads 

W [GHZn) := ^1 - \GHZ N )(GHZ N \, (19) 

which has also been used in an experiment This 
witness detects genuine iV-qubit entanglement around an 
7V-qubit GHZ state. The constant 1/2 corresponds to the 
maximal squared overlap between the GHZ state and the 
pure biseparable states, this can be calculated using the 
methods presented in Ref. [l2| . 

Witness Eq. Q19[l can be interpreted in the following 
way. (W ( - GHZn) ) is -1/2 only for the GHZ state state. 
For any other state it is larger. In general, the more 
negative {W {GHZn) )\^) is, the closer state |*) is to the 
GHZ state. It is known that in the proximity of the GHZ 
state there are only states with genuine ./V-qubit entan- 
glement, so the constant in Eq. (|19|l is chosen such that 
if (W<- GHZn) ) < then the state is in this neighborhood 
and is detected as entangled. 

Now we will show that the projector-based witness 
Eq. (|19fl can be constructed as a sum of elements of the 



stabilizer. First, note that the generators of the stabi- 
lizer g^ HZN ^ can j-jg use d to define a very convenient 
basis for calculations, namely their common eigenvec- 
tors. \GHZm) is one of them giving a +1 eigenvalue 

for all Sj^ HZ s. Allowing both +1 and a —1 eigenval- 
ues in Eq. 2^ ./V-qubit states can be defined which 
are orthogonal to each other and form a complete ba- 
sis. These are all GHZ states, but in different bases, i.e., 
they are of the form {\x (1 \..x {N ^)± \x [1) ...x {N) ))/y/2 with 
G {0, 1} and x^ = 1 - x^ . We will refer to this ba- 
sis as the GHZ state basis. This basis will turn out to be 
extremely useful since all operators of the stabilizer are 
diagonal in this basis. 

Now, the projector onto a GHZ state can be written 
with the stabilizing operators 0] 

N MGHZn) ^ 

\GHZ N ){GHZ N \ = \\— (20) 

k=l 

This can be directly seen in the GHZ basis. From 
Eq. H2U|) it follows that the projector-based witness can 
also be decomposed into the sum of stabilizing terms, 
i.e., Y\}( GHZn ) is a stabilizer witness. This decomposition 
can be used to measure the witness locally. However, the 
number of settings needed seems to increase exponen- 
tially with the number of qubits [l5]. Thus we have to 
find other stabilizer witnesses for which measuring them 
is more feasible. 

For detecting biseparable entanglement it was enough 
to use two stabilizing operators for the witness. For the 
detection of iV-qubit entanglement we have to make mea- 
surements on all qubits and have to measure a full set of 
generators. This is because we need that the expectation 
value of the witness is minimal only for the GHZ state. 
If the witness does not contain a full set of generators 
then there are at least two of the elements of the GHZ 
basis giving the minimum. There is, however, always a 
superposition of two basis vectors of the GHZ basis which 
is biseparable. This biseparable state would also give a 
minimum for our witness. Thus this witness could not 
be used for detecting genuine iV-qubit entanglement. 

The main idea of detecting genuine multi-qubit entan- 
glement with the stabilizing operators is the following: 
Observation 2. Let us consider some of the stabilizing 
operators. If this set of operators contains a complete set 
of generators and for a given state the expectation val- 
ues of these correlation operators are close enough to the 
values for a GHZ state, then this state must be close to 
a GHZ state and is therefore multi-qubit entangled. 

Now we can derive the first entanglement witness. 
Theorem 2. The witness 

N 

W:=(N-1)1-J2sI GHZn) (21) 

k=l 

detects only states with genuine N-qubit entanglement. 
Proof. In Eq. (|21|l the constant term, cq = N — 1, was 
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chosen such that the observable 



# of qubits 



X a := W - aW (GffZjv) > 0. 



(22) 



for some a > is positive semidefinite. Then we have for 
any state g the inequality aTr{gW {GHZN) ) < Tr(gW) 
which implies that Tr(gW) < can only happen if 
Tr{gW^ GHZN) < 0. Thus all states detected with W are 
also detected by W( GHZ «) and W is a multi-qubit wit- 
ness. Clearly, we would like to have Co as small as possi- 
ble, since the smaller cq is, the more entangled states W 
detects. Simple calculation leads to cq = N — 1. One can 
check that with thischoice W - 2W^ GHZn '> > 0. □ 
Both witnesses W (GffZjv) and W detect entangled 
states close to GHZ states. The main advantage of W 
in comparison with Y\}( GHZn > lies in the fact that for im- 
plementing it only two measurement settings are needed, 
namely the ones shown in Fig. ^a). The first setting, 

{IW.IP),..,!^)}, is needed to measure (s[ GHZn) ), 
the second one, is to measure the 

expectation values for the other generators. The other 
characteristic to check is the noise tolerance of the wit- 
ness. The witness W detects states mixed with noise of 
the form Eq. I|14|) if p ao ise < 1/N. Thus the noise toler- 
ance decreases as the number of qubits increases. 

However, using a similar construction it is also possible 
to derive a witness which is robust against noise even 
for many qubits and still requires only two measurement 
settings: 

Theorem 3. The following entanglement witness detects 
genuine N-qubit entanglement for states close to an N- 
qubit GHZ state 



yy{GHZ N ) ._ 31 _ 2 



(GHZ N ) 



1 



^k 



n 

k=2 



1 



This witness has the best noise tolerance among stabilizer 
witnesses which need only two measurement settings and 
have the property yy( GHZ «) - 01 W { - GHZn) > for some 
a>0. 

Proof. To prove the first statement, one can show by di- 
rect calculation that W {GHZn) - 2W ( - GHZn) > 0. Thus 
yy(GHZ N ) j g & mu lti-qubit witness. For the proof of op- 
timality please see Appendix B. □ 
To give an example, for the simple case of three qubits 
the witness is 

3 



W 



(GHZ 3 



■l-lWlt 2 '!' 3 ' 



I[ Z (D Z (2) +Z (2) Z (3) +Z (1) Z (3)]_ 



(24) 

Threc-qubit genuine multi-qubit entangled states can be- 
long to the so-called W-class or to the GHZ-class 16J. 
Knowing that >V (GHZ3) + 1/4 detects GHZ-class entan- 
glement we obtain that yv( GffZ3 ) + 1/2 detects also 
only GHZ-class entanglement. 



c/3 
tofj 

g 

I 

<^ 
O 



2mmmmmm 

(a) 

1 [x] [z] [x] [z] [x] [z 



2 \z\ [x] \z\ [x] \z\ ^x 

(b) 



FIG. 1: (a) Measurement settings needed for the approach 
presented in this paper for detecting entangled states close 
to GHZ states and (b) cluster states. For each qubit the 
measured spin component is indicated. 



The witness Y\?( GHZn ) is quite robust against noise. It 
detects states mixed with white noise as true multipartite 
entangled for 



< 



3_2( 2 -W)' 



(25) 



Thus it tolerates at least 33% noise, independently from 
the number of qubits. Again, only two measurement set- 
tings are necessary for an implementation [see Fig.^a)]. 

The expression for the witness Y\}( GHZn ) can be sim- 
plified using the fact that 



V g(GHZ N ) 



n 

k=2 



1 



= |00...0) <00...0| + |11...1)(11...1|, 



(26) 

where |11...1)(11...1| is the projector on the state with 
all spins down and |00...0)(00...0| is the projector on the 
state with all spins up. Using this one obtains 

w (ghz n ) _ 21 _ X (2) X (3) ■ ■ ■ A (w) 

-2|00...0)(00...0| - 2|11...1)(11...1|(27) 

E. Bell inequalities for GHZ states 

As a sidestep we will discuss now a very surprising 
feature of the stabilizer witnesses, namely that they are 
closely related to Mermin-type Bell inequalities [jj], [3 • 
As we will see, this relationship sheds new l ight on the 
question whether and when Bell inequalities 120L l2l) 
can detect genuine multipartite entanglement. 

First note that witnesses different from the previous 
ones can be obtained by including further terms of the 
stabilizer and using more than two measurement settings. 
For instance, following the lines of the previous section 
it is easy to see that the observable 



21 - s[ GHZs) [t + S^^p + S 3 



(GHZ 3 ) 



21 + yWy( 2 )x (3) + x (1) y (2) y (3) 

y(i)x( 2 )y( 3 ) - A (1) A (2) A (3) 



(28) 
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detects genuine three-party entanglement around a GHZ 
state. It detects a GHZ state mixed with white noise 
if Pnoisc < 1/2. The witness WMermin is equivalent to 
Mermin's inequality [Tsj 

(iWi( 2 'i (3) ) - (x (1) r (2) r (3) ) - 

(y(l) X (2) F (3) ) _ ( y(Dy(2) x (3) ) < 3 (2g) 

In Ref. the condition given in Eq. 12911 was used for 
detecting entanglement in a 3-qubit photonic system and 
a measurement result equivalent to (WMermin) ~ — 0.83± 
0.09 was obtained, thus the created state was genuine 
three-qubit entangled. 

For N > 3 the Bell operator in Mermin's inequality 
contains also only stabilizing terms: 
Theorem 4. For the Bell operator of the Mermin's in- 
equality fU 

1 

y(l)y(2) X (3) X (4) 



M N := 



y(l)y(2)y(3)y(4) . , . x (N-l) X (N) _ ^ 



(30) 

the operator expectation value for biseparable states is 
bounded by 



(Mn) < -, 



(31) 



while the quantum maximum is 1. Note that a term in 
Eq. I.Vfl)) represents the sum of all its possible permuta- 
tions. 

Proof. Mn can alternatively written as 

N g(GHZ N ) 



M N = S[ 



(GHZ* 



n 

k=2 



1 



= |00...0)(11...1| + |11...1)(00...0|. (32) 

The maximum for (Mn) for biseparable states can be 
obtained knowing 

\GHZ N )(GHZ N \ - M N = \GHZ N ){GHZ N \ > 0, (33) 



where \GHZ 
biseparable states g 



(|0000...> - |1111...))/V2. Hence for 
1 



(M N ) e < (\GHZ N )(GHZ N \) e < 



For full y se parable states the bound is lower and is given 
in Ref. 



(34) 

en 

□ 



III. WITNESSES FOR CLUSTER, GRAPH AND 
W STATES 

A. Witnesses for cluster states 

Let us turn now to cluster states. These are a family 
of multi-qubit states which have attracted increasing at- 
tention in the last years. A cluster state |CV) is defined 



to be the state fulfilling the equations \Cn) 
with the following stabilizing operators 



sI Cn) \c n ) 



q(Cn 



z (k-l) x (k) z (kA 



;fc = 2,3,...,JV-l, 



(35) 



"■JV 

Witnesses similar to Eq. (|I1|I which rule out full 
separability can be constructed with two locally non- 
commuting operators as pfj 



W (C N ) ._ t _ S (CN)_ s (C* 



i ^ - /, for fc- 1,2, ...,7V- 1. (36) 

This witness detects biseparable states close to an N- 
qubit cluster state. The proof is essentially the same as 
the one for Eq. lfTT|l . 

Note that Wi involves only two generators which 
act on at most four-qubits. This witness detects whether 
the reduced density matrix of the qubit quadruplet cor- 
responding to qubits (fc — 1), (k), (fc + 1) and (k + 2) is 
entangled. The state of the rest of the qubits does not 



influence (W ( k CN> ). The witness W k 

7-moise ^ 1/2- 

The following witnesses have a better noise tolerance 



,(Cn) 



tolerates noise if 



W, 



(Cm) 



1 



q 



for k = 1,2,. 



fc+i 
N- 1 



S. 



(Civ) q{Cn) 



(37) 



This witness still involves only the qubits of a quadruplet 
and tolerates noise if p n oisc < 2/3. 



Using witnesses W 



(CV) 



one can construct a "compos- 



ite" entanglement witness for which the noise tolerance 
increases with the number of qubits: 
Theorem 5. The following entanglement witness detects 
entangled states close to a cluster state 



K-l 



yw(Cjv) . = "IT i _ $( Cn 

' v comp ' 11 



4fc+l 



_ c(Cn) _ c(Cjv) q(C N ) /OCA 
°4fc+2 °4fc+l°4fc+2' \°°) 



k=0 



where K := [(-/V + 2)/4j and \_x\ denotes the integer part 
of x. The witness Eq. ]3^l tolerates noise if 



Pu 



< 



2 K + r 



(39) 



Proof. The noise tolerance comes from direct calculation 
along the lines of Eq. (|15|1 . Note that all the terms in 
the product in Eq. Ij38(l act on disjoint sets of qubits. 
This is the reason that such a composite witness can be 
constructed. □ 
One can also construct witnesses for the detection of 
genuine multipartite entanglement close to cluster states, 
similar as for the case of GHZ states: 
Theorem 6. The witnesses 



yw(CV) ._ 

2 



31 



Cn){Cn 

n 



(40) 



q(Cn) 



+ 



n 



s 



(Cn) 



odd k 



(41) 



detect genuine N-party entanglement close to a clus- 
ter state. yy( Clv ) is optimal from the point of view of 
noise tolerance among the stabilizer witnesses which need 
only two measurement settings and have the property 
yy(Civ) _ a >v(Civ) > o f or some a > 0. 
Proof. First we have to prove that Eq. (|40|l is a multi- 
qubit witness. For that we have to use that from a cluster 
state one can generate a singlet between arbitrary qubits 
by local operations [|| . Using the results of Ref. |27| this 
implies that the maximal Schmidt coefficient of a clus- 
ter state when making a Schmidt decomposition with re- 
spect to an arbitrary bipartite split does not exceed the 
maximal Schmidt coefficient of the sin glet , which equals 
l/y/2. Then, from the methods of Ref. [IJ it follows that 
yy(cv) j g a w it ness f or multi-qubit entanglement. After 
that we have to prove that Eq. (|41 1) is a multi-qubit wit- 
ness. This can be proved similarly as it has been done for 
Theorem 3 using that W {Cn) - 2W {Cn) > 0. Concerning 
the optimality, see Appendix B. □ 
The stabilizing operators in W^™) are again di- 
vided into two groups corresponding to the two settings 
{X^,Z^,X^,Z^\...}and {Z^\X^ 2 \Z^\X^\...} 
as shown in Fig. ^b). Witness W^ Cn ' tolerates mixing 
with noise if 

J (4-4/2TT)- 1 for even N 

Pnoisc <{ [4 _ 2(1/2^ + 1/2^)]-! for odd N ' 

(42) 

Thus, for any number of qubits at least 25% noise is 
tolerated. 

B. Witnesses for graph states 

Results similar to the ones derived before hold also 
for graph states 0, . These states are defined in the 
following way: One takes a graph, i.e., a set of N vertices 
and some edges connecting them. Edges of this graph are 
described by the adjacency matrix T. Tki = 1 (0) if the 
vertices k and I are connected (not connected) . Based on 
that one can define the stabilizing operators 

Sfrt-xW [] Z^> =X^Y[(Z^f*'. (43) 

Neighbors Z^fc 

I of k 

Then, the graph state \Gn) is defined as the 7V-qubit 
state fulfilling \G N ) = S [Gn) \G n ). Physically, the graph 
provides also a possible generation method: The graph 
state can be created from a fully separable state by us- 
ing Ising interactions between the connected qubits. In 
fact, many useful multi-qubit states can be treated in the 
graph state formalism, for instance also GHZ states and 
cluster states. The corresponding graphs are shown in 
Fig.0 

Theorem 7. A witness detecting biseparable entangle- 
ment close to graph states can be given as 

wif N) := 1 - S { h GN) - SI GN \ (44) 
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1 2 




(a) (b) (c) 



FIG. 2: Graphs corresponding to different graph states, (a) 
Star graph and (b) cluster state graph. A graph state corre- 
sponding to a star graph is equivalent to a GHZ state under 
local unitaries. (c) A seven- vertex graph which has a triangle. 
Due to this it is not a two-colorable graph. 

where the spins (k) and (I) are neighbors, and a witness 
detecting genuine N-party entanglement can be defined 
as 

w (g n ) : =(N-l)t-J2 S k GN) - ( 45 ) 
k 

Proof. The proofs are essentially the same as before. 
First one has to show that 1/2 — \Gn)(Gn\ is a witness 
for true multipartite entanglement. Then one can prove 
that witnesses Eq. I|44|l and Eq. (|45|l detect also only 
genuine multipartite entanglement. □ 
Now the number of settings needed for measuring 
yy(G N ) depends on the graph corresponding to the graph 
state to be detected. For this we need the notion of 
colorability of graphs. A graph is M-colorable, if one 
can divide the vertices into M groups and assign to the 
vertices of each group a color, such that neighboring ver- 
tices have different colors (see Fig. 2). For two-colorable 

graphs, only two settings are needed. In this case the 
(c ~t 

>Sfc operators can be divided into two groups corre- 
sponding to the two settings. In general, for M-colorable 
graphs M settings are needed for measuring the witness 
yy(Giv)_ j n ^his senses, the most settings (M = N) are 
needed for the complete graph |29l |. 

It is, however, important to note at this point that the 
colorability of the graph is not an intrinsic and physical 
property of the graph state. Usually, a graph state can 
be represented by different graphs up to local unitaries, 
i.e., different graphs G\ and G2 can result in two graph 
states \G\) and \G%) which are the same up to a local 
change of the basis. Here, G\ and G2 may have differ- 
ent colorability properties, e.g., G\ may be two-colorable 
and G2 iV-colorable. The question, which graphs give the 
same graph state is still an open and challenging problem 
in stabilizer theory. Recently, much progress has been 
achieved concerning a subclass of local unitary transfor- 
mations, the so-called local Clifford transformations [30| . 
Efficient algorithms have been developed which allow the 
determination of all graphs which are equivalent within 
this subclass of local unitary transformations. These 
methods can readily be used to find witnesses with a 
small number of measurement settings. 
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C. Obtaining the fidelity of the prepared state 

Let us say that in an experiment we intend to create 
a GHZ state. Beside knowing that the prepared state p 
is entangled, we would also like to know how good its 
fidelity is. The fidelity could be measured by measuring 
the projector on the GHZ state 



F :=Tr(\GHZ N )(GHZ N \p). 



(46) 



However, we would encounter the same problem as with 
witnesses: The number of local settings needed for mea- 
suring the projector increases rapidly with the size. 

Fortunately it is possible to obtain a lower bound on 
the fidelity from the expectation value of our witnesses. 
For example, for our GHZ witness defined in Theorem 3 
we have yv(GHZ N ) _ 2 y^{GHZ N ) > g Hence 

\GHZ N )(GHZ N \ > i - 1 W {GHZn \ (47) 

where W( GHZ ») is defined in Eq. 1(23)1. Now a lower 
bound on the fidelity F := Tr{\GH Z N )(GH Z N \p) can 
be obtained as F' := Tr(P'p) where 



p> ~ I _ \yj(GHZ N ) 
' 2 2 

Q (GHZ N ) fl N Q (GHZr. 

_ £l +_± , TT 2k 

9 11 9 



- 1. (48) 



k=1 



Note that for measuring P' only two local measurement 
settings are needed. 

Let us see how good this lower bound is for the noisy 
GHZ state g(p n oisc) defined in Eq. 1(15(1 . For this state we 
have 

F (Pnoisc ) 1 .Pnoisc (1 2 ) , 

F'( Paoisc ) = l-p noisc (3/2-2-( Ar - 1 )). (49) 

The difference is AF w j? n oise/2 for large N. Bounds can 
be obtained similarly for the fidelity with respect to the 
cluster state based on 

\C N ){C N \>\-\w^\ (50) 



D. Witnesses for mixed states 

In an experiment, after detecting entanglement in the 
prepared state, one might be interested to measure and 
drop one of the qubits and investigate the state obtained 
this way. Here we assume that we do not know the mea- 
surement result thus the system is in a mixed state corre- 
sponding to the reduced density matrix of the remaining 
qubits. We will show with an example that witnesses can 
easily be derived also for this mixed state. 



Let us consider a concrete example, namely, measuring 
a qubit of the four-qubit cluster state: 
Example 1. The witness 



W ( 83 ) . = 1 _ z <i) z <?) _ xWxWzW. 



(51) 



detects entanglement around the state g 3 which is ob- 
tained from the four-qubit cluster state 



\C' A ) := |0000) + |0011) + |1100) - |1111) 



(52) 



after the fourth qubit is measured. Note that for \C'±) we 
used a local transformation in order to be able to present 
the cluster state in a convenient form in the Z basis. 
Proof. After measuring the fourth qubit, the three-qubit 
mixed state is 



Q3 := 

id := 

ID := 



1 

2 

1 

71 
i 

V2 



(ID(D + ID(D), 

(|00) + |ll))<8)|0) ) 

(loo) - lii» ® li). 



(53) 



Note that g 3 is biseparable, i.e., the third qubit is unen- 
tangled from the first two. 

Now we will determine the stabilizer of g 3 . The sta- 
bilizers of states l^), given with their generators, are 
:= (ZWzt^lWlf^^.T^)). We would like 
to detect entangled states close to a particular mixture of 
these two states. Any mixture of states |^) is stabilized 
by the common elements of 5^ ) : The stabilizer of these 
mixed states is = (ZWzW,XWxWzW). 

Based on that, a stabilizer witness Eq. (|51|l can now be 
constructed. The constant in Eq. (|51|l was determined 
such that W^ 3 ) detects indeed entangled states only. 
This can be proved similarly to the derivation of The- 
orem 1. Witness W^ e ^ needs two measurement settings 
and tolerates noise if p n oise < 1/2. □ 



E. Entanglement witnesses with nonlocal 
stabilizing operators 

There are states which do not fit the stabilizer for- 
malism, however, it is still possible to find a simpler 
witness than the one obtained by decomposing the pro- 
jector. As an example let us look at the W state 
\W 3 ) = (| 100) + |010) + |001»/V3. For this state, the 
projector-based witness is known to be 



W im) = -1 
3 



\W 3 )(W 3 



(54) 



It tolerates noise if p no ise < 8/21 ss 0.38 and measuring 
a local decomposition of this witness requires five mea- 
surement settings In this section we will present 
witnesses for the \W 3 ) state which need only three and 
two measurement settings, respectively. 
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What is new with the stabilizing operators of IW3)? 
Clearly, we have to leave the requirement that these op- 
erators should be a tensor product of single-qubit oper- 
ators. Now the stabilizing operators must be allowed to 
be the sum of several such locally measurable terms. 

A set of stabilizing operators with simple local decom- 
position can be found in the following way. Let us assume 
that we create the \Wz) state from state |000) using a 
unitary dynamics U . The generators of the stabilizer for 
1 000) are 



g(|000)) ._ z (k) 



for k = 1,2,3. 



(55) 



Hence one can get the generators for a group of operators 
which stabilize \Wa) 



s (w 3 ) = us \ooo) u i > 



(56) 



Let us try to find U. It must fulfill \W 3 ) = U\000), i.e., 
written in the product basis we must have 



1 

7! 



[uiu 2 ...us\ 



Hence, u\ can be obtained. The other columns of U 
are not determined by Eq. (|57J) and constrained only by 
requiring that U is unitary (i.e., the columns of the ma- 
trix must be orthonormal to each other). Thus U is not 
unique. A possible choice for an U fulfilling Eq. 1571) is 

U := -±-(X^Z^ + X^Z^ + Z^X^). (58) 
V3 

The generators of a group of stabilizing operators can 
be obtained based on Eq. (|5oTl 



(w 3 



~,(W 3 . 



( Z (D + 2Y^Y^Z^ + 2X^Z^X^), 



:= -{ZM + 2Z^Y^Y^ + 2X (1 ^X^Z^), 
3 

:= ^(Z& + 2Y^Z^Y^ + 2Z^X^X {i) ). 

(59) 



These three stabilizing operators uniquely define the W 
state. Again, by multiplying them with each other, other 
stabilizing operators can be found. However, with the 
exception of s[ m) S^ m) = , they are 

all nonlocal. 

Now let us try to create an entanglement witness which 
detects genuine multi-party entanglement, but is still eas- 
ier to measure than the projector-based witness W^ 3 ). 



Consider the following witness constructed from some of 
the elements of the group generated by S 



(w 3 ) 



(W 3 ) Q (W 3 ) 



w = Co t-sr 3 's: 



j(w 3 ) s (_w 3 ) 



s (w 3 ) s (w 3 ) 



'2 -3 , (60) 

where Co and c\ are positive constants. The expectation 
value of W is minimal for the W state. Constants cq and 
Ci must be determined such that for some a > we have 
W — aW^ 3 ' > 0, and also we have the possible best 
noise tolerance. Thus we have: 

Theorem 8. A witness detecting genuine three- qubit en- 
tanglement around a [W^) -state is 
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w (w 3 



-1 + 2Z^Z^Z^ 



2X (k) x (i) 



2 y( fe )yW_zWzW).(6i) 



{k)y(l)\ 



It requires three measurement settings and tolerates noise 

if Pnoise 

< 4/15 « 0.27. 
One can also leave out the Z^ terms: The following 
entanglement witness detects genuine three-qubit entan- 
glement close to a IW3) -state 



( 5? ) W':=(l + V5)l 



X^X (2) - A (2) A (3) - A (1) A (3) 

y(l)y(2) _ y(2)y(3) _ y(l)y(3)_ 

(62) 



This witness requires the measurement of two settings 
and tolerates noise if p n oisc < (3 — VS}/4 w 0.19. The 
proof is given in Appendix C. 



IV. CRITERIA USING VARIANCES AND 
UNCERTAINTY RELATIONS 

Let us now explain, how criteria in terms of variances 
can be derived, using uncertainty relations 0, |U • 
In Ref. the following recipe was presented for the 
derivation of such criteria, called the local uncertainty 
relations (LURs). Consider a bipartite quantum system 
and let Ai be some observables on one party, fulfilling a 
bound 



(63) 



for all states on this party. Here, S 2 (Ai) = (A 2 ) — (A) 2 
denotes the variance of the state. This inequality is an 
uncertainty relation for the Ai, with Ua > iff the ob- 
servables Ai have no common eigenstates. Let us assume 
that we have also observables Bi on the second system, 
fulfilling a similar bound 5 2 {Bi) > Ub- Then, we may 
look at the observables Mi = Ai £g> 1 + 1 £g> Bi on the com- 
posite system. As it was shown in Ref. [33 for separable 
states 



J2S 2 (M t )>U A + U B 



(64) 
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has to hold, and a violation of this bound implies that 
the state is entangled. 

Now we show how the witness Wm HZN ^ in Eq. 1111) 
for GHZ states can be improved using nonlinear terms: 
Theorem 9. Let us define 

A x = XWXW.-.XM, 
A 2 = Z (k \ 

Bl = _ I (w) X M,.. I m i 

B 2 = -Z^ k+1 \ (65) 

fork= 1, 2, N— 1. Using these operators, the following 
necessary condition for separability can be given 

i - {sf HZN) ) - (s£Jf Zn) ) 

^((A 1 +B 1 ) 2 + {A 2 +B 2 ) 2 )>0. (66) 

Proof. For the uncertainties of observables Ak and Bk, 
one has the bounds 

S 2 (A 1 ) + S 2 (A 2 ) > 1, 

(5 2 ( J B 1 )+^ 2 ( J B 2 ) > 1. (67) 

Knowing that (X^) 2 + (Z (fe) ) 2 < 1, which implies that 
5 2 {X^) + <5 2 (Z( fe >) > 1 these bounds should not be a 
surprise. A detailed proof which relies on the fact that 
A\ and A 2 anti-commute, is given in Appendix C. From 
this and the method of the LURs Eq. I|66|) follows. □ 
Eq. (|66[) is a nonlinear necessary condition for separa- 
bility. It can be considered as a nonlinear "refinement" 
of W^ HZn) in Eq. (TTJ, which improves the witness 
by subtracting the squares of some expectation values. 
The fact that LURs can sometimes improve entangle- 
ment witnesses was, for bipartite systems, first observed 
in Ref. |34j . There, the magnitude of the improvement 
for a special case was also investigated numerically. In 
the given case for GHZ states it is important to note that 
the LUR does not improve the noise tolerance, when the 
GHZ state is mixed with white noise. This is due to the 
fact that for the totally mixed state as well as for the 
GHZ state the squared mean values in Eq. I|66|) vanish. 
However, there are many states which are not of this 
form and which are detected by the LUR and not by the 
witness. 

In fact, many of the witnesses from the previous sec- 
tions can be improved via the method presented above. 
For instance, one may add the extra observables A% — 

x^xw ■ ■ ■ xV'-VyW , b 3 = -r( fe + 1 )x< fc + 2 ) • • • it") 

to the observables A\, ...,B 2 from above. Then S 2 (Ai) + 
5 2 {A 2 ) + S 2 (A 3 ) > 2 and the same bound holds for the 
Bi. This leads to the separability criterion 

- ^{{A 1 +B 1 ) 2 + {A 2 +B 2 ) 2 + {A- i + B 3 ) 2 ) >0 

(68) 



which improves the witness of Eq. (|13fl . Also the witness 
in Eq. (|36|) can be improved using the same methods, 
leading to the separability condition 

1 - {S^ N) )-{S^ + f)-\{Z^X^ -Z^) 2 - 
_I( Z (fc)_x( fe+1 >z( fe + 2 )) 2 >0, (69) 

and the witness for graph states can be improved in a 
similar manner. 

To conclude, it turned out that LURs can improve the 
presented witnesses which were designed for ruling out 
full separability. The witnesses for the detection of gen- 
uine multipartite entanglement can not be so simply im- 
proved, mainly for two reasons. First, LURs are specially 
designed for bipartite systems, and no generalization to 
multipartite systems in known so far. Second, accord- 
ing to our definition of multipartite entanglement a state 
which is a convex combination of states which are bisep- 
arable with respect to different partitions is also bisep- 
arable. This implies that it is not enough to show that 
a state is biseparable with respect to each partition in 
order to show that it is multipartite entangled. How- 
ever, if one defines a state to be multipartite entangled 
if it is not biseparable with respect to any partition (as 
it is done sometimes in the literature, see, e.g., Ref. |35j) 
then LURs can be used to detect multipartite entangle- 
ment by ruling out biseparability for every bipartite split 
(as proposed in |32|). 



V. CRITERIA USING ENTROPIC 
UNCERTAINTY RELATIONS 

Let us now discuss another possibility of deriving sepa- 
rability criteria in terms of uncertainty relations, namely 
criteria based on entropic uncertainty relations. As we 
will see, the stabilizer formalism allows us to formulate 
easily such criteria, however, they are not as strong as 
the witnesses or the variance based criteria. 

The recipe to derive such criteria was described in 
Ref. (3j| and goes as follows: If we have an observable 
M we can define V(M) e — (p x , ...,p n ) as the probability 
distribution of the different outcomes for measuring M 
in the state g. One can characterize the uncertainty of 
this measurement by taking the entropy of this probabil- 
ity distribution, i.e., by defining H{M) := H(P(M) e ). 
Here, we only consider the entropy to be the standard 
Shannon entropy H(V) := — J2kPk^ n (Pk), however, a 
more general entropy like the Tsallis entropy may also 
be used. Given two observables M and N which do not 
share a common eigenstate, there must exist a strictly 
positive constant C such that H(M) + H(N) > C holds. 
The difficulty in this so-called entropic uncertainty rela- 
tion (EURs) lies in the determination of C. For results 
on this problem see Refs. [3(| |37j ■ For the detection of 
entanglement, the following result has been proved [33| : 
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Let A\ , A 2 and B\ , B 2 be observables with nonzero eigen- 
values on Alice's (respectively, Bob's) space obeying an 
EUR of the type 

H(Ai) + H(A 2 ) > C (70) 

and the same bound for B\,B 2 . If g is separable, then 

H{A 1 ®B 1 ) e + H{A 2 ®B 2 ) e >C. (71) 

holds. For entangled states this bound can be violated, 
since A\ ® B\ and A 2 <8> B 2 might be degenerate and have 
a common entangled eigenstate. 

In order to apply this scheme to the detection of en- 
tanglement in the stabilizer formalism we have to recall 
some more facts. Assume that we have two observables 
of the form M = Mi Pi an d N = V\ v jQj where the 
Pi,Qj are the projectors on the eigenspaces. Here, we do 
not require M and N to be non-degenerate, i.e., the P{ 
and Qj may have ranks larger than one. In this situation, 
it was shown in Ref. [37j that for these observables the 
EUR 

H(M) + H(N) > - In max (72) 

ij 

holds, where = \\PiQj\\ = 

max|^) yj (ip\(PiQj)^ (PiQj)\i(j) is the norm of the 
operator PiQj. This has two consequences. First, it 
follows immediately that for one qubit 

H(X)+H(Y)>\n(2), (73) 

holds, and similar relations hold for Y, Z etc. Second, 
if A is an observable consisting of Pauli measurements 
on N qubits, e.g., A = X^Y^...Z^ N \ then for N + 1 
qubits the EUR 

H(A <g> X {N+1) ) + H(1 2 n <g> Y {N+1) ) > ln(2), (74) 

holds. Here 1 2 n denotes the identity on the first N 
qubits. Similar result hold also, if the observables on 
the qubit N + 1 are replaced by other Pauli matrices. 
Eq. (|74ll can directly be proved from Eq. I|72[l by identi- 
fying the corresponding Pi and Qj . 

Armed with these insights, we can formulate now en- 
tropic criteria for stabilizer states: 

Theorem 10. For GHZ and cluster states, respectively, 
the following necessary conditions for biseparability can 
be given using entropic uncertainties 

N 

^H(Si GHZ ^) > ln(2), (75) 

N 

Y, H ( S l CN) ) > M 2 )- ( 76 ) 

fc=i 

Any state violating these conditions is genuine N-qubit 
entangled. Note that for GHZ and cluster states the left 
hand side of Eq. j75| ) and Eq. j76] ), respectively, is zero. 



Proof. To prove Eq. i|75|) it suffices to look at pure 
biseparable states, since the entropy is concave in the 
state. So let us assume that \if}) = \a)\b) is a biseparable 
state. For definiteness, we assume that \a) is a state of 
the the qubits 1,2, ...,k and \b) is a state of the qubits 
k + 1, k + 2, N; the proof of the other cases is similar. 
In order to apply the recipe from above, we have to show 
that on the first k qubits the EUR 

k 

H{X^...X^) + Y J H{S i h GHZN) ) + H(ZW) > ln(2) 

i=2 

(77) 

holds. This is easy to see, since H(X^ ■ ■ ■ X^) + 
H(Z (k y) > ln(2) is valid, due to Eq. d). Similar ideas 
can be used for cluster states. The proof is essentially 
the same as for the GHZ state. Similar statements can 
also be derived for arbitrary graph states. □ 

As we have seen, it is quite straightforward to for- 
mulate entropic criteria for stabilizer states. However, 
one should clearly state that the criteria in the presented 
form are not very strong. For instance, the criterion in 
Eq. (|75|l detects states mixed with white noise only for 
Pnoiso < 0.123 for N = 3 and p noisc < 0.083 for N = 4. 
The robustness to noise can, as shown for some two-qubit 
cases in Ref. |33|. be improved, if other entropies than 
the Shannon entropy are used. However, for these en- 
tropies no such general bound as in Eq. (|72|l is known. 
So only a better understanding of the EURs can help to 
explore the full power of the entropic criteria. 



VI. CONCLUSIONS 

In summary, we have shown that stabilizer theory can 
be used very efficiently to derive sufficient criteria for 
entanglement. Knowing some stabilizing operators of a 
state allows for an easy derivation of a plethora of entan- 
glement criteria which detect states in the vicinity of the 
state. This holds for linear as well as nonlinear criteria, 
and for the different types of entanglement in the multi- 
partite setting. We also noted that the resulting criteria 
exhibit several interesting features: They all are easy to 
implement in experiments, some of them have interest- 
ing connections to Bell inequalities, others are nonlinear 
improvements of witnesses. 

A natural continuation of the present work lies in the 
systematic extension of the presented ideas to states 
which do not fit directly in the stabilizer formalism. 
GHZ, cluster and graph states are not the only multi- 
partite states which are interesting from the viewpoint of 
quantum information science. As we have shown in the 
example of the W state, also for states outside the stabi- 
lizer formalism similar ideas can be applied by identifying 
their nonlocal stabilizing operators. Exploring this direc- 
tion might help to clarify the structure of multipartite 
entanglement. 
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APPENDIX A: The stabilizer group 

We summarize the properties of the stabilizing oper- 
ators of a given A-qubit quantum state \^f). These are 
locally measurable operators S for which S\$f) = 
They have ±1 eigenvalues. For states considered in this 
paper (GHZ, cluster and graph states) they are tensor 
products of Pauli spin matrices. 

Let us now list the properties of the set S containing 
the stabilizing operators of For any Si, S2, S G S 

SiS 2 e S, 
[Si, S2] = 0, 

S 2 = 1. (78) 

Now it is clear that the elements of S form a commutative 
(Abelian) group. It is called the stabilizer. One of its 2 N 
elements is given as 



N 



s 



k ■- 



Y[(SiT kl 



(79) 



i=i 



where au is the 2th digit of the binary number (i.e., 
iV-tuple of {0, 1}) corresponding to the number k £ 
{0, 1, 2, 2 N — 1}. Operators {S;}^ are the generators 
of the group. If ctki = 1 then we say that Sk contains 
generator Si in its definition of the type Eq. H79|) . 



APPENDIX B: Proof of optimality for the GHZ and 
cluster state witnesses 



In this section we will prove that the GHZ and clus- 
ter state witnesses defined in Eq. (|23*|) and Eq. 
respectively, are optimal. That is, it is not possible to 
find a stabilizer witness W which needs only two mea- 
surement settings, has the property that for some a > 
we have W — aW > and has better noise tolerance than 
the witnesses presented in this paper. Here W denotes a 
projector-based witness. 



Before presenting the proof, let us analyze what 
we understand on measurement settings. Let us con- 
sider operators which can be measured with one mea- 
surement setting. These form a group, members of 
which commute locally (for an explanation of local com- 
muting see Sec. II. B). The projector cannot be mea- 
sured with two settings, since the stabilizing opera- 
tors cannot be divided into two locally commuting sub- 
groups. So we will find two subgroups, which con- 
tain as many stabilizing operators as possible. Let 
us consider GHZ states first. The largest such sub- 



(GHZ N ) Q (GHZ N ) 
' a 3 



■' N 



Here the 



group is (S 

group is given by the generators. The operators in 
this group commute locally since they all contain only 
Z( fe )'s, and neither X( fc )'s nor FW 's. All the other ele- 
ments of the stabilizer have the form s[ GHZn ^Q where 
Q G {sf HZN \sf HZN \...,S {GHZN) ). All these opera- 
tors contain only X^'s and Y^'s, and do not contain 
J^( fe )'s. None of these commute locally with any other 
element of the stabilizer, except with the identity. Thus 
the other subgroup can contain only one such operator 
and the identity. 

Now the two subgroups corresponding to the two local 
settings are 



£^GHZm) 
^•{GHZm) 



I cj{GHZn 



/q(GHZn) q(GHZn) 

\°2 5 °3 J 



,s 



(GHZpt)\ 
N I 



(80) 



Any operator which is a linear combination of operators 
in £1 can be measured by the first measurement setting. 
Similarly, any operator which is linear combination of 
operators in £2 can be measured by the second measure- 
ment setting. It is an important question whether there 
is another choice for C\ and £2- In order to answer this, 
we have to remember that we want to detect genuine 
A-qubit entanglement. Thus we must find a witness for 
which (W) has a unique minimum for GHZ states. For 
that, we must be able to measure a complete set of gen- 
erators with the two settings. (For more details, see Sec. 
II. D. ) It can be seen that no other two subgroups can 
be found which fulfill this requirement. For simplicity, 
in the following we choose Q — 1 . The two measurement 
settings corresponding to this case are shown in Fig.^a). 
Choosing Q to be not the identity would not change our 
argument and would not lead to witnesses with better 
noise tolerance than the ones presented here. 

As discussed before, the eigenvectors of the genera- 
tors {sI^ HZn ^}^ =1 of the stabilizer form a complete ba- 
sis (GHZ basis) . We will use it to represent states of the 
A-qubit Hubert space. Let us use A-tuples of {0, 1} for 
labeling the basis states. If the A;th digit is (1) then 



for the basis state ^g^ HZN ^ 



-1 ((Sf HZnI ) = -!)• 



Now I [00...0]) is the GHZ state. Here square brackets are 
used in order to draw our attention to the fact that the 
GHZ basis is not a product basis and the N digits do not 
correspond to a physical partitioning of the system. 
Next, let us use this labeling to order the 2 N ba- 
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sis states from | [00...0]) to |[ll...l]). What is the ma- 
trix form of 5"( GffZjv ) m this basis? It is clearly di- 
agonal. Moreover, it must give +1 and —1 expec- 
tation values for states of the form |[0s2S3---Sjv]) arid 
|[1s2S3...sjv]), respectively. Knowing this, it must have 
the form diag(+l, —1) (8 l 2 iv-i. Here diag denotes a di- 
agonal matrix and the size of the identity is indicated 
by a subscript. Note again that the tensor product does 
not correspond to a physical partitioning of the system. 
However, the matrix form of s[ GHZn ' in the GHZ basis is 
the same as the matrix form of Z^ in the product basis. 
Similarly, the matrix form of g£* HZN > ^ e QHZ basis 
is the same as the matrix form of Z^ fe ' in the product 
basis. 

After these considerations about measurement set- 
tings, let us start our proof. Let us find out how op- 
erators measurable by the first setting look like in the 
GHZ basis. Now it is clear that such operators must 
have the form A £g> where matrix A is of size d\ = 2 
and e?2 = 2( N ~ 1 \ Operators measurable by the second 
setting have the form 1^ eg) B in this basis. Here B is 
of size d 2 ■ We construct our witness from two operators 
corresponding to the two measurement settings as 



W := cl, 



A®1 



1« 



B, 



(81) 



where c is a constant. Now we will find an operator W for 
which for some a > we have W — aVV^ GHZN ' > > and 
which is optimal from the point of view of noise tolerance. 
Without the loss of generality, we set a — 2. Both W 
and the projector witness y\?( GHZN ) are diagonal in the 
GHZ basis. Hence from the condition W > 2W {GHZn) , 
the following constraints on the diagonal elements and 
bk of A and B, respectively, can be obtained 



c — ai — bi > —1, 

c-dk-k > +1; k > 2; I > 1, 

c-ak-k > +1; k > 1; I > 2. 



(82) 



The maximal noise pn m it tolerated by our witness can be 
computed as given in Eq. I|15|) . For this formula we have 
to use 



{GHZ N \W\GHZ N 

-N 



c — a\ — b\ 



(83) 

In order to find the optimal witness, constant c and 
the elements of A and B must be chosen such that 
the noise tolerance pumit is maximized, under the con- 
straints Eq. (|82|l . This is the case if all the three inequal- 
ities in Eq. ijS^ are saturated. Thus we obtain c = 3, 
a\ — b\ — 2, dk = bk = for k > 2 and 



1 



Plimit 



2/di - 2/d 2 



(84) 



The witness obtained this way coincides with the witness 
given in Eq. (jSSJ- 

Now let us turn to cluster states. We will need the 
following lemma: 

Lemma 1. Let us consider the following three subsets of 
the stabilizer S^ Cn '> 



Vy := {A e S {Cn) : A contains S- t 



, c(C N ) 



(85) 



(C \ 

and does not contain S\ , x }, 



V 2 := {A e S {Cn) : A contains Sfc* 



(C«) 



(86) 



and does not contain S\ N }, 
V 3 := {A e S {Cn) : A contains both 
and S^ N) }. 

(87) 

If a local measurement setting makes it possible to mea- 
sure an operator in V n , then it does not make it possible 
to measure any operator in V m for m ^ n. 
Proof. First let us prove our Lemma for I = 1. Table QI a) 
shows for some elements of the stabilizer which measure- 
ments on qubits (1) and (2) are needed to measure them. 
Here for brevity superscript Cat is omitted. A binary 
pattern indicates whether a given element of the stabi- 
lizer contain or does not contain generators S%, S2, and 
53. For example, entry 100 represents operators which 
contain Si and do not contain S2/3. Such operators are, 
for example, Si, S1S4 and SiS^Sq. In the second row of 
the table "A","F" and " Z" represent Pauli spin matri- 
ces. " 1" indicates that no measurement is needed on the 
given qubit. For example, in the column of 101 entry 
"XI" indicates that for measuring SiS 3 = X^X^Z^ 
an X measurement is needed for qubit (1) and no mea- 
surements are needed for qubit (2). 

The left two columns correspond to operators of set Vi 
which contain Si and do not contain 5*2 . The middle two 
columns correspond to operators of set Vi which contain 
S2 and do not contain S\. The right two columns corre- 
spond to operators of set V3 which contain both Si and 
S2. Based on the table we see that only 100 and 101 have 
compatible measurements for the first two qubits: "XZ" 
and "XV, respectively. An operator of V n cannot be 
measured together with any operator of V m for m =/= n. 
Thus Lemma 1 works for 1 = 1. 

Now let us prove, that if our Lemma is true for I = 
k — 1, then it is true for I — k. Table Hfb) shows the 
measurements for qubits k and k + 1 for a particular 
set of operators. An entry of the table, let us say 0100, 
represents all operators which contain Sk, and does not 
contain Sk-i,Sk+i and Sk+2- Now the measurements for 
qubits k and fc+1 for 0101 seem to be compatible with the 
measurements for 1010, 1011, 1110, and 1111. However, 
we assumed that our Lemma is true for I = k — 1. Thus 
0101 and 1010 cannot be measured with the same setting. 
For the same reason, 0101 and 1011 cannot be measured 
with the same setting. Taking the other operators having 
compatible measurements for qubits k and fc + 1, we find 
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S1S2S3 


100 


101 


010 


Oil 


110 


111 


(D (2) 


xz 


XI 


zx 


ZY 


YY 


yx 



(a) 



Sk-iSkSk+iSk+2 


0100 


0101 


1100 


1101 


Q(k)Q{k + l) 


xz 


XI 


YZ 


Yl 


Sk-iSkSk+iSk+2 


0010 


0011 


1010 


1011 


Q(k) (k + 1) 


zx 


ZY 


IX 


1Y 


Sk-iSkSk+iSk+2 


0110 


0111 


1110 


1111 


Q(k)Q{k + l) 


YY 


YX 


XY 


XX 



(b) 



TABLE I: (a) Local measurements needed for measuring some 
elements of the stabilizer on qubits (1) and (2) for a cluster 
state. The binary pattern indicates presence or absence of 
generators Si, S2 and 53. (b) Local measurements needed for 
qubits k and k + 1 where k > 2. The binary pattern indicates 
presence or absence of generators Sk~i, Sk, Sk+i and Sk+2- 
For details see text. 

that measuring them with the same setting is prohibited 
by our Lemma for I = k — 1. □ 
Lemma 1 limits the size of the subgroup of the sta- 
bilizer of \Cn) measurable by one measurement setting. 
Namely, for even N the largest subgroup can be given 
with N/2 generators, while for odd N the maximum is 
(N + l)/2. Two subgroups of the stabilizer allowing for 
the measurement of a complete set of generators are 

/*(Cn) / c(Cn) q(Ojv) q(Cn) \ 

L -l ■— Wl j J 3 ! D 5 > •■■II 

£< c «) ._ ^g(cv) g{c N ) gCcv) ^ (88) 

Now di=d 2 = 2 N ' 2 for even N and d 1/2 = 2 iN±1 ^ 2 for 
odd N. There are other possibilities for the two measure- 
ment settings, however, they do not give a better noise 
tolerance since the optimal noise tolerance Eq. (|84J) de- 
pends only on d\ and d 2 , and does not depend on which 
elements of the stabilizer are measured by the two set- 
tings. Optimization leads to witness (|41|l constructed 
using two projectors, corresponding to the two settings. 
Thus witness given in Eq. (|41|l is also optimal. 

APPENDIX C: Some technical calculations 

Calculations for the W state — Here, we prove that 
Eq. (|62[1 describes an entanglement witness detecting 
states with three-qubit entanglement around a W state. 
Note that this witness was constructed independently 
from the stabilizer theory. The proof is also useful in 
general, since it shows how to find the minimum of an 
operator expectation value for biseparable states analyt- 
ically. 

Let us first assume (1)(23) biseparability. Then for a 
state of the form * = ® * 2 3 

(W') = (1 + y/E) - - (X«) 

-(x&>xW) - (yU)(yW) - <yW)<yW) 

_/y(2)y(3)\ 

= (1 + V5)-(X^) 1 [(X^+X^) 23 ] 

+ ( y(i) )l[( y(2) + y(3) )23] + [(X (2) X (3) + Y ^Y^) 23 ] 



where 

z:=(* (1) )i; y-=(Y^h. (90) 

Here and (...)23j respectively, denote expectation 

value computed for ^1 and ^23- To be explicit, matrix 
F, as the function of two real parameters, is given 

F(x,y) := (l + V5)l-x[X (2) +X (3) ]-?;[Y (2) +Y (3) ] 
- [XWX& (91) 

Note that operator F acts on the second and third qubits 
while its two parameters depend on the state of the first 
qubit. The expectation value of F(x, y) can be bounded 
from below 

(F(x,y)) 23 > A min [F(x,y)} =y/E- yj\ + 4(.x 2 + y 2 ), 

(92) 

where A m ; n (F) is the smallest eigenvalue of F. The right 
hand side of Eq. (|92|) was obtained by finding the matrix 
eigenvalues analytically. Using (XW) 2 + (Y~W) 2 < 1, 
we obtain that (W) is bounded from below by 

(W) > min {A min [F(x,y)}} = 0. (93) 

x +3/ <1 

Thus our witness has a non-negative expectation value 
for biseparable pure states with partitioning (1)(23). Due 
to that W' is symmetric under the permutation of qubits, 
(W) > holds also for pure biseparable states with a 
partitioning (12) (3) and (13) (2). It is easy to see that 
the bound is also valid for mixed biseparable states. In 
contrast, for three-qubit entangled states we have (W) > 
V5 — 3. Among pure states, the equality holds only for 
\W 3 ) and IW3) = (|011) + |101) + |110))/V3, and their 
superpositions. 

Calculations for the LURs described in Sec. IV — 
To compute the bounds required for the derivation of 
the LURs we show the following. Let Aj, for i — 
1, 2, n be some observables, which anti-commute pair- 
wise (i.e., AiAj + AjAi = for all i 7^ j) and which 
have all ±1 eigenvalues (i.e., Af — l for all i). Then 
2™=i ( ^ 2 ( j 4i) ^ n — 1. To show this, is suffices to show 
that Y^i=i(Ai) 2 < 1- This can be proved as follows: 
We take real coefficients Ai,...,A n with J2™ =1 \ 2 = 1. 
Then, using the fact that the Ai anti-commute, we have 

(Eti a^) 2 = ELi x i A i = H=i = a - So ' fOT a11 

states KE"=i x i A i)\ = IEr=i^(A)| < 1 holds. Since 
the Aj are arbitrary, this implies that the vector of the 
mean values ((^-1), (A 2 ),..., (A n )) has a length smaller 
than 1. Thus, YJi=i\ A i) 2 < 1 has to hold - Tb - is method 
can be used to derive all the bounds in Sec. IV. 
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